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Abstract 

Measurements of the muon charge asymmetry in inclusive pp —> WX production at ^fs = 1 TeV are presented. 
The data sample corresponds to an integrated luminosity of 4.7 fb -1 recorded with the CMS detector at the LHC. With 
a sample of more than twenty million W —> pv events, the statistical precision is greatly improved in comparison to 
previous measurements. These new results provide additional constraints on the parton distribution functions of the 
proton in the range of the Bjorken scaling variable x from 10 -3 to 10 _1 . These measurements are used together with 
the cross sections for inclusive deep inelastic ep scattering at HERA in a next-to-leading-order QCD analysis. 


1. Introduction 

In the proton-proton collisions at the Large Hadron 
Collider (LHC), the primary processes for inclusive 
W-boson production are the annihilation processes: 
ud —> W + and du —> W ~, where a valence quark comes 
from one proton and a sea quark comes from the other 
proton. Due to the presence of two valence u quarks 
as opposed to only one valence d quark in the proton, 
more W + bosons are crated compared to W~ bosons. 
The measurement of the asymmetry in the production 
of W + and W~ bosons as a function of boson rapidity 
can provide constraints on the d/u ratio and the sea 
quark densities in the proton. 

There are large uncertainties on the parton distri¬ 
bution functions (PDFs) for the scale and Bjorken x 
region explored in the p-p collisions at LHC at ^[s = 1 
TeV corresponding to Bjorken x from 10~ 3 to 10 _1 . In 
addition, there are differences between different PDF 


Email address: saranya. ghosh@cern. ch (Saranya Samik 
Ghosh (for the CMS Collaboration)) 

1 Speaker 


sets. A high precision measurement of the W-boson 
charge asymmetry can contribute to the improvement 
of the knowledge of PDFs. 

It is difficult to measure the W-boson production 
asymmetry in terms of the boson rapidity because of 
the presence of neutrinos in the leptonic W decay. 
Instead, the lepton charge asymmetry in terms of 
lepton pseudorapidity is measured. The lepton charge 
asymmetry is defined as 

¥(W + -> tv) - ¥(W~ -> tv) 

sm = £ - p- - (i) 

%(W + -> £ + v) + %(W- -» e-v) 
where is the differential cross section for W- 

dr] 

boson production and subsequent leptonic decay and 
rj = -ln[tan(6/2 )] is the charged lepton pseudorapidity 
in the laboratory frame, with 6 being the polar angle 
measured with respect to the beam axis. 

The measurement of lepton charge asymmetry in 
the muon channel, corresponding to the W —> pv decay 
channel of the W-boson, using a data sample corre¬ 
sponds to an integrated luminosity of 4.7 fb -1 recorded 
with the Compact Muon Solenoid (CMS) detector 0 
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is presented along with a next-to-leading order QCD 
analysis using the CMS measurement combined with 
the cross sections for deep inelastic e + ~p scattering at 
HERA (3). 

2. Measurement of muon charge asymmetry 

The events selected are required to contain one well 
reconstructed isolated muon with transverse momentum 
( pr ) greater than 25 GeV. An additional measurement 
of muon charge asymmetry is also made with a higher 
muon p T threshold of 35 GeV. Events containing a 
second isolated muon with p T greater than 15 GeV are 
rejected, in order to reduce the Drell-Yan background. 
The important backgrounds constitute of Drell-Yan 
events, QCD multijet events with muons in the final 
state, W —> tv with the r decaying into a muon and 
neutrino and tt events. Monte Carlo (MC) simulated 
samples, that include a full simulation of the CMS 
detector, are used to help evaluate the background 
contributions in the data sample. 
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Figure 1: The dimuon mass profile as a function of muon r/ for p+ 
after the correction. The reference is obtained using MC generator 
level information and taking into account the effect of reconstruction 
on the resolution. 


shown in Figure [2] 


The selected events are divided into eleven bins 
based on the absolute pseudorapidity (\r]\) of the muon, 
corresponding to [0.0-0.2], [0.2-0.4], [0.4-0.6], [0.6- 
0.8], [0.8-1.0], [1.0-1.2], [1.2-1.4], [1.4-1.6], [1.6-1.85], 
[1.85-2.1], [2.1-2.4] and the muon charge asymmetry is 
calculated separately in these eleven bins of muon \tj\. 

Corrections are applied in order match the pileup, 
that is multiple proton-proton interactions in the same 
event, in MC to that in data (the MC simulation is 
generated with a different pileup distribution than that 
we observe in the data) and also to match the boson 
transverse momentum distribution in MC to that in 
data. The muons momentum scale is corrected for 
any bias due to misalignment and mismodelling of the 
magnetic field. Corrections are based on the charge, 
pseudorapidity and the azimuthal angle of the muons 
and have been derived on Z —> pp selected sample. The 
effect of muon momentum scale corrections on positive 
muons is shown in Figure [l] 

To account for the Drell-Yan background, the 
Drell-Yan MC sample is normalized with theoretical 
cross-section to luminosity of the data sample. MC to 
data scale factors are used to account for the difference 
in efficiency of muon identification and selection 
between data and MC. Additional k-factors have been 
calculated to correct normalization. The correction 
factors are obtained using Z —> pp samples in bins of 
dimuon invariant mass. The Drell-Yan normalisation is 


x10 3 CMS, L = 4.7 fb' 1 at Vs = 7 TeV 



- 2-10 12 

ri of leading muon 


Figure 2: The 77 distribution of the leading p+ in Z/y* —> p + p 
sample. The dimuon invariant mass is within 60 < m m <120 GeV. 
The MC simulation is normalized to the data. 

The missing transverse energy (MET) distribution 
is asymmetric in 0 due to detector and reconstruction 
effects. Corrections are applied to correct for MET 0 
asymmetry. Even after correcting for the 0 asymmetry, 
the distribution of MET in MC is different from that in 
data. MET is used for signal extraction, and the MET in 
MC is matched to that in data using recoil corrections. 
Recoil corrections are corrections to the hadronic recoil, 
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which is the vector sum of transverse momenta of all 
particle candidates excluding the candidate muon. The 
corrections are derived on a Drell-Yan selected control 
sample. The MET distribution in the Drell-Yan selected 
sample in data and MC after applying the corrections 
is shown in Figure [3] The corrections are verified in a 
QCD control sample that is obtained by inverting the 
muon isolation criteria. The MET distribution in the 
QCD control sample is shown in Figure [4] 


xIO 3 CMS, L = 4.7 fb' 1 at\/s = 7 TeV 
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Figure 3: Data to simulation comparison for MET in the DrellYan 
control sample. 
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Figure 4: Data to simulation comparison for MET in the QCD control 
sample for p+ events. The W —> pv contribution (light shaded region) 
is normalized to the integrated luminosity of the data sample using a 
MC simulation, and the normalization of the QCD simulation (dark 
shaded region) is taken as the difference between the data and the 
estimated W —> pv contribution. 


timation is done using the MET distribution in data and 
MC samples. MET for W + and W~ candidate events are 
taken separately for the different \r/\ bins used and MET 
templates with corrections from MC fitted to MET dis¬ 
tribution in data using binned maximum likelihood fit. 
Simultaneous but separate fits for W + and W~ candi¬ 
date events is performed in each \rj\ bin with the yields 
of W + —> yu + v, W~ —> /Tv and QCD normalization 
floated in the fits. A total of 12.9 million W + and 9.1 
million W candidate events give yields of ~ 84% for 
W —> pv, ~ 8 % QCD, ~ 8 % electro weak background 
consisting of Drell-Yan and W —> tv and a small frac¬ 
tion of tt. An example of the MET fits performed for 
the muon p T > 25 GeV sample for pseudorapidity bin 
0.0 < I 77 I < 0.2 is shown in Figure [5] 
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Figure 5: Fits using MET distributions for the muon pj > 25 GeV 
sample for the extraction of the W —> pv signal from data for pseu¬ 
dorapidity bin 0.0 < \r]\ < 0.2. The fits to W + —> p + v (a : top) 
W~ —» p~v (b : bottom) are shown separately. 


The process of signal extraction and background es- The asymmetry obtained after the signal extraction 
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process is then corrected for difference in efficiency of 
selection of positive and negative muons. 


3. Results 

The systematic uncertainties associated with the 
measured muon charge asymmetry are calculated. The 
largest source of systematic errors comes from the 
errors on the measured efficiencies of p + and jT , that 
are used to correct the asymmetry for difference in 
efficiency of p + and p~. Another large source of sys¬ 
tematic uncertainty comes from the QCD background 
estimation. The error from the QCD background 
estimation is taken from two parts, the error coming 
from the ratio of positive and negative muon events 
in the QCD background and the uncertainty in the 
shape of the MET distribution coming from the QCD 
background. Some of the other major sources of sys¬ 
tematic uncertainty come from uncertainty in the muon 
momentum scale, effect of final state radiation (FSR) 
(FSR is not corrected for in the measured asymmetry 
but systematic uncertainty is assigned to account for 
its effect) and PDF uncertainties coming from the MC 
samples used in the measurement. Figure [6] displays 
the uncertainties coming from different sources for the 
Pt > 25 GeV measurement of muon charge asymmetry. 



Figure 6 : Square of the uncertainties from different sources added in 
quadrature for uncertainties associated with the pr > 25 GeV mea¬ 
surement of muon charge asymmetry for the different \rj\ bins. 

As a cross-check, muon charge asymmetry is mea¬ 
sured separately for the positive and negative rj regions 


by performing an identical measurement in 22 muon r) 
bins. The asymmetry obtained from the negative and 
the positive regions is compared as shown in Figure [ 7 ] 
for the p T > 25 GeV case. The charge asymmetries for 
tj > 0 and p < 0 regions are in good agreement with 
each other. 



Figure 7: Comparison of the muon charge asymmetry extracted 
for the positive pseudorapidity (77 > 0 ) and negative pseudorapidity 
(77 < 0) regions for muon pj > 25 GeV sample. The uncertainties 
include only the statistical uncertainty from the signal extraction and 
uncertainty in the determination of the efficiencies for positive and 
negative muons. 


The final measurement of the muon charge asym¬ 
metry along with the uncertainties is compared to 
predictions by five different PDF sets as shown in 
Figure [5] for the muon p T > 25 GeV sample and in 
Figure [9]for the muon p T > 35 GeV sample. The PDF 
predictions are calculated at NFO using NFO FEWZ 
3.1. The measured muon charge asymmetry is in good 
agreement with predictions from CT10, NNPDF and 
HERA PDF sets. The agreement with predictions from 
MSTW PDF set is poor, although it is slightly improved 
for MSTW2008CPdeut PDF set. 


A comparison of the measured muon charge asym¬ 
metry result to the previous CMS electron charge 
asymmetry measurement 0| extracted using 0.84 fb -1 
of the 2011 CMS data is shown in Figure [lO] The 


lepton p T threshold on the electron channel measure¬ 
ment is 35 GeV and the comparison is made with the 
corresponding muon channel result. The muon and the 
electron charge asymmetry are consistent with each 
other. 
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Figure 8: Comparison of the measured muon charge asymmetry 
for the muon pj > 25 GeV sample to the NLO predictions calcu¬ 
lated using the FEWZ 3.1 MC tool interfaced with the NLO CT10, 
NNPDF2.3, HERAPDF1.5, MSTW2008, and MSTW2008CPdeut 
PDF sets. 


CMS, L = 4.7 fb' 1 at\fs = 7 TeV 



Figure 9: Comparison of the measured muon charge asymmetry 
for the muon pj > 35 GeV sample to the NLO predictions calcu¬ 
lated using the FEWZ 3.1 MC tool interfaced with the NLO CT10, 
NNPDF2.3, HERAPDF1.5, MSTW2008, and MSTW2008CPdeut 
PDF sets. 


CMS, \[s = 7 TeV 
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Figure 10: Comparison of measurement of muon charge asymmetry 
to the previous CMS electron charge asymmetry result. Results are 
shown for lepton pj > 35 GeV. 


4. QCD PDF Analysis 

Combination of CMS muon charge asymmetry mea¬ 
surement and deep inelastic scattering (DIS) data from 
HERA is used do determine the PDFs of the proton 
in an NLO perturbative QCD (pQCD) analysis. The 
procedure for the determination of the PDFs follows 
the approach used in the HERAPDF1.0 QCD fit (3). 
The PDF uncertainties take into account experimental, 
model, and parametrization uncertainties according to 
the general approach of HERAPDF1.0. 

A comparison of the PDFs obtained for valence u 
and d quarks from fits using HERA data and muon 
asymmetry measurements to that using HERA data 
only is shown in Figure m The muon charge asym¬ 
metry measurements, together with HERA DIS cross 
section data, improve the precision of the PDF of 
valence quarks over the entire v range. 

5. Conclusion 

The W —> fiv lepton charge asymmetry is measured 
in proton-proton collisions at yfs = 7 TeV using a data 
sample corresponding to an integrated luminosity of 
4.7 fb -1 collected with the CMS detector at the LHC. 
The asymmetry is measured in 11 bins in absolute 
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Figure 11: Distributions of u valence (top) and d valence (bottom) 
quarks as functions of x at the scale Q 2 = . The PDFs obtained 

from fit to the HERA data and muon asymmetry measurements (light 
shaded band), and to HERA only (dark hatched band) are compared. 
The total PDF uncertainties are shown. In the bottom panels the dis¬ 
tributions are normalized to one for a direct comparison of the uncer¬ 
tainties. The change of the PDFs with respect to the HERA-only fit is 
represented by a solid line. 


muon pseudorapidity, \tj\, for two different muon p T 
thresholds of 25 and 35 GeV. 

The measurement has significantly lower statisti¬ 
cal and systematic uncertainties compared to previous 
CMS measurements. The total uncertainty per bin 
is 0.2 to 0.4%. The asymmetry measures in data is 
in good agreement with the theoretical predictions 
using CT10, NNPDF2.3, and HERAPDF1.5 PDF 
sets. The agreement with the prediction based on the 
MSTW2008 PDF set is poor, although the agreement is 
improved when using the MSTW2008CPdeut PDF set. 
With smaller uncertainties compared to current PDF 
uncertainties, this measurement can be used to signif¬ 
icantly improve the determination of PDFs in future fits. 

An NLO QCD analysis using the muon charge 
asymmetry along with DIS data from HERA improves 
the prediction of the valence quark distributions from 
range x = 10 -4 to 0.5. 
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